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技术 方法 


基于 sift 特 征 点 的 透射 电子 显微镜 多 幅 图 像 拼接 
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摘要 :目的 为 了 解决 TEM 在 高 分 辩 率 下 仅 观 察 到 小 视野 的 问题 ,本 文 提 出 了 一 种 基于 sift 特 征 点 的 TEM 多 幅 图 像 拼接 方法 。 


方法 首先 ,提取 待 拼接 图 像 的 sift 特 征 点 。 然 后 ,利用 sift 特 生 


F 点 结合 多 项 式 拟 合 校正 场 的 方法 对 图 像 进行 几何 校正 ,再 利用 


sift 特 征 点 进行 图 像 配 准 。 最 后 ,采用 泊 松 图 像 编 辑 的 方法 进行 融合 以 达到 无 颖 拼接 。 结 果 60 张 TEM 局 部 肾 小 球 图 像 能 拼接 
出 完整 的 肾 小 球 ,图 像 的 配 准 误差 为 62.5~187.5 nm, 且 与 其 它 3 种 图 像 拼 接 方法 相 比 ,能 够 有 效 减 小 形变 .避免 伪 影 的 产生 。 结 
论 本 文 提 出 的 方法 能 得 到 高 分 辩 率 、 大 视 时 的 TEM 图 像 ,为 肾 活 检 病 理 诊断 提供 有 力 的 支持 。 

关键 词 :图 像 拼 接 ; 特 征 点 ;图 像 几何 校正 ;图 像 配 准 ; 泊 松 图 像 编 辑 ; 高 分 辨 率 ; 大 视野 


Multiple transmission electron microscopic image stitching based on sift features 


LTMz LU Yanmeng, HAN Shuaihu’, WU Zhuobin’, CHEN Jiajing',, LIU Zhexing’,, CAO Lei’ 
‘School of Biomedical Engineering, Electron Microscope Lab, School of Basic Medical Sciences, Southern Medical University, Guangzhou 


510515, China 


Abstract: We proposed a new stitching method based on sift features to obtain an enlarged view of transmission electron 


microscopic (TEM) images with a high resolution. The sift features were extracted from the images, which were then combined 
with fitted polynomial correction field to correct the images, followed by image alignment based on the sift features. The 


image seams at the junction were finally removed by Poisson image editing to achieve seamless stitching, which was validated 
on 60 local glomerular TEM images with an image alignment error of 62.5 to 187.5 nm. Compared with 3 other stitching 
methods, the proposed method could effectively reduce image deformation and avoid artifacts to facilitate renal biopsy 


pathological diagnosis. 
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透射 电子 显微镜 (TEM) 检 查 是 凤 活 检 病 理 检 查 的 
重要 组 成 部 分 , 肾 组 织 内 伴 特 殊 有 形 结构 的 各 种 肾脏 疾 
病 的 诊断 及 其 新 疾病 的 发 现 ,大 多 通过 电镜 检查 得 以 确 
诊 "。 北 京 大 学 肾脏 病 研究 所 对 777 例 肾 活 检 标 本 的 病 
理 诊断 结果 表明 ,电镜 对 18.5% 的 病例 起 决定 性 诊断 作 
用 ,对 13.5% 的 病例 诊断 提供 重要 参考 作用 ”。 但 是 用 
于 透射 电镜 检测 的 样品 要 求 标本 块 很 小 ,通常 超 注 切片 
中 只 含有 1~2 个 肾 小 球 ,本 身 对 于 反映 病变 程度 已 有 一 
定 的 局 限 性 ,再 者 目前 透射 电镜 普遍 配备 的 高 分 状 率 
CCD 数字 化 相机 的 视野 都 比较 小 ,一 张 图 片 只 能 包含 
病变 的 局 部 ,缺乏 反映 广泛 病变 区 的 信息 ,容易 造成 漏 
诊 。 如 果 要 包含 完整 的 病变 区 ,如 一 个 完整 的 凤 小 球 ， 
就 只 能 在 较 低 的 放大 倍率 下 进行 图 像 采 集 , 才 失 了 电子 
显微镜 高 分 辩 率 的 优势 ,拍摄 的 图 像 分 辨 率 难 以 满足 临 
床 诊断 的 需要 。 这 就 产生 了 一 个 矛盾 : 既 要 保证 图 像 清 
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晰 又 要 包含 完整 的 病变 区 。 目 前 凤 活 检 所 采集 的 TEM 
图 像 很 难 达 到 两 全 ,大 大 限制 了 TEM 在 临床 超 微 病理 
活检 中 的 应 用 ”。 人 解决 此 问题 最 直接 的 方法 就 是 将 多 
张 高 分 辩 率 且 有 一 定 重 车 区 域 的 TEM 图 像 进 行 配 准 后 
拼接 出 一 幅 高 分 辨 率 、 大 视野 的 图 像 ,为 临床 医生 提供 
完整 病变 区 域 的 高 分 辨 率 几 像 以 提高 病变 诊断 效率 。 

图 像 配 准 的 方法 主要 分 为 基于 像素 的 配 准 和 基于 
特征 的 配 准 。 基 于 像素 配 准 的 方法 是 利用 两 幅 图 像 重 
全 区 域 的 灰 度 信息 构建 相似 性 度量 函数 ,求解 匹配 参 
数 ,从 而 建立 图 像 的 配 准 。 常 用 的 基于 像素 的 配 准 方法 
有 :基于 相关 性 ` 互 信息 以 及 频 域 傅 里 时 方法 等 “。 该 类 
配 准 方法 能 充分 利用 图 像 信 息 ,但 对 重合 区 域 小 的 图 像 
很 难 佑 计 准 确 的 匹配 参数 ”。 基 于 特征 的 配 准 方法 涉 
及 到 点 特征 \ 线 特征 和 区 域 特征 *, 点 特征 是 最 基本 的 特 
征 , 其 适应 范围 最 广 生 对 其 计算 及 描述 都 相对 简单 “， 
此 ,本 文选 择 特征 点 的 配 准 。 基 于 特征 点 的 方法 通过 提 
取 具 有 几何 不 变性 的 点 ,由 特征 点 之 间 的 匹配 得 出 匹配 
参数 ,从 而 建立 图 像 的 配 准 。 几 何不 变性 的 特征 点 及 其 
描述 子 的 获取 是 该 类 方法 的 关键 。Smith 等 "用 
SUSAN 算 子 提取 角 点 ,但 对 角 点 的 定位 不 够 准确 。 
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Harris 等 提取 出 了 具有 高 重复 性 检测 率 的 harris 角 
点 。Schmid 等 "在 此 基础 上 ,加 入 方向 不 变性 的 描述 
子 ,使 该 角 点 具有 方向 不 变性 的 特点 ,但 并 不 具备 尺度 
不 变性 。Lowe 等 "在 高 斯 差分 空间 提取 出 了 对 图 像 的 
尺度 .方向 具有 不 变性 是 有 一 定 仿 射 不 变性 的 sift 特 征 
点 。Yan 等 中 对 sift 特 征 点 的 描述 子 用 PCA 方 法 进行 降 
维 得 到 pca-sift 特 征 点 ,虽然 减少 了 计算 量 ,但 与 sift 特 
征 点 相 比 ,在 尺度 变化 大 以 及 图 像 较 模 糊 时 ,匹配 效果 
较 差 。Mikolajczyk 等 将 sift 特 征 点 中 描述 子 转 换 为 
对 数 极 坐 标 中 计算 ,得 出 gloh 特 征 点 ,虽然 与 sift 特 征 点 
相 比 ,gloh 特 征 点 在 光照 变化 大 及 图 像 较 模糊 时 ,匹配 
更 有 优势, 但 计算 量 更 大 "”。Bay 等 “利用 hessian 算 阵 
来 提取 surf 特 征 点 ,提高 了 运算 速度 且 在 图 像 较 模 糊 时 
优 于 sift 特 征 点 ,但 所 估计 的 匹配 参数 不 够 稳定 。Juan、 
Kayning 等 将 这 三 种 特征 点 进行 比较 ,指出 sift 特征 
点 在 尺度 .旋转 不 变性 以 及 匹配 纹理 复杂 图 像 中 均 
优 于 以 上 三 种 特征 点 ,在 仿 射 不 变性 中 sift 与 surf、gloh 
相当 。 

由 于 sift 特 征 点 具有 较 好 的 尺度 和 旋转 不 变性 ,对 
复杂 纹理 图 像 的 配 准 ,所 估计 的 参数 也 相对 稳定 ,因此 
本 文选 择 基于 sift 特 征 点 的 配 准 方法 对 纹理 复杂 的 肾 小 
球 TEM 图 像 进行 配 准 。 由 于 显微镜 成 像 系统 导致 了 图 
像 的 形变 和 亮度 不 均匀 ,故而 本 文 分 别 进行 了 配 准 的 预 
处 理 和 后 处 理 。 本 文 第 1 节 介 绍 了 st 特征 点 的 提取 、 
两 幅 TEM 图 像 的 拼接 过 程 以 及 多 幅 图 像 的 完整 拼接 过 


程 ,第 2 节 是 实验 结果 ,第 3 节 为 本 文 讨论 。 


1 方法 
1.1 sift 特 征 点 的 提取 上 

sift 特 征 点 是 指 图像 中 具有 尺度 方向 不 变性 上 且 有 
一 定 仿 射 不 变性 的 点 ,提取 方法 如 下 : 

第 1 步 : 构 建 图 像 高 斯 空间 和 高 斯 差分 DOG 
(Difference of Gaussian) 空 间 , 如 图 1A 所 示 。 

令 I(x,y) 为 输入 图 像 , 由 式 (1) 得 到 高 斯 空间 第 一 
组 第 一 层 图 像 Li(x,y, o0): 

Li(x,y,00)= G(x,y, oo0FI(x, y) (1) 

G(x,y,00) 表 示 尺 度 为 0 的 高 斯 函数 。* 污 示 卷 积 。 

然后 由 公式 (2) 构 建 高 斯 空间 的 第 "组 第 ; 层 图 像 
L(x;y, 0): 

L(x,y,0%)= Gx, y, oFL, (x,y, 0 )， 
(0'=ko’, ',o= JJ -(0 ) ,i=2,3,..., m,n=1,...,1,) 
其 中 ,1 为 高 斯 空间 总 组 数 ,m 为 每 组 图 像 总 层 数 ， 
人 = 2 儿 - 为 尺度 系数 。 

再 由 公式 (3) 构 建 DOG 空间 中 第 zx 组 第 ; 层 图 像 
D' (x, y， 0): 


(2) 


Dy) Cpe Ty h 

(=1,2,...,m-)) 

本 文中 取 o, =1.6,k= V2,m=5,1=4。 

第 2 步 :提取 特征 点 ,并 确定 特征 点 方向 。 特 征 点 
在 DOG 空 间 中 同 组 是 相 邻 的 三 层 图 像 中 提取 ,如 图 1B 
所 示 , 若 e 所 示 像 素 点 的 灰 度 值 是 周围 26 像 素 点 灰 度 值 
的 极 值 , 则 将 该 像素 点 作为 特征 点 。 若 c 为 特征 点 所 在 
图 像 层 的 尺度 ,以 特征 点 为 中 心 ,计算 3 半径 范围 内 像 
素 点 的 梯度 ,如 图 1C, 再 对 该 范围 内 梯度 进行 直方 图 统 
计 , 将 统计 出 的 主 方向 作为 特征 点 的 方向 ,如 图 1D。 

第 3 步 :计算 特征 点 的 128 维 向 量 描述 子 。 本 文 在 

特征 点 周围 16 x 16 令 域内 计算 描述 子 。 先 计算 出 该 区 
域内 梯度 的 方向 和 幅 值 ,为 保持 方向 不 变性 ,将 该 区 域 
内 所 有 像素 的 方向 以 特征 点 的 方向 为 基准 进行 调整 。 
然后 将 该 邻 域 分 解 为 4 x 4 个 子 区 域 。 最 后 将 各 子 区 域 
中 的 梯度 统计 到 8 个 方向 中 ,获得 128 维 向 量 描述 子 。 
1.2 图 像 拼 接 
图 像 拼 接 是 利用 多 张 有 一 定 重 到 范围 的 小 视野 的 
图 像 ,获取 一 幅 高 分 辩 率 无 颖 的 大 视野 图 像 的 方法 。 
本 文 图 像 拼 接 的 方法 由 图 像 几何 畸变 校正 、 图 像 配 准 、 
图 像 融合 三 部 分 组 成 ,其 中 基于 sift 特 征 点 的 图 像 配 准 
是 图 像 拼接 的 关键 。 
1.2.1 图 像 配 准 "" 本 文 基于 sift 特 征 的 图 像 配 准 方法 
包括 以 下 3 步 :第 1 步 :sift 特 征 点 的 匹配 。 令 图 像 i、j 
为 相 邻 的 两 张 图 像 , 根 据 1.1 节 的 方法 分 别提 取 图 像 
i 、j 中 sift 特征 点 集 X ,={xY ,Xx2, .XX = 
{Xxx .… x}。 对 图 像 中 的 某 特征 点 x0 的 描述 子 ， 
在 X 的 描述 子 中 用 BBF(Best Bin First) 算 法 搜索 与 
其 欧式 距离 最 近 和 次 近 的 sf 特征 点 x 和 x 中 ,并 计算 
出 最 近 距 离 和 次 近 距 离 的 比值 。 当 该 比值 小 于 0.8 时 ， 
将 x ,x 加 入 匹配 点 集 。 最 终 得 到 图 像 ;中 的 特征 点 集 
i ee ,其 在 图 像 7 中 对 应 点 集 为 
六 ={2 x 站, 且 (x ,x ) 为 一 组 匹配 寺 征 点 ,n 为 
特征 点 总 匹配 数 ;第 2 步 :图 像 匹 配 。 假 定 图 像 i、j 满 足 
如 下 匹配 模型 : 


(3) 


0) @ my hy @ 
Xx = Hx = hy hx (4) 
0 1 


其 中 ,HH 是 一 个 仿 射 变换 矩阵 ,可 根据 RANSAC 
(Random Sample Consensus) 算 法 计算 出 ;第 3 步 :将 图 
像 投影 到 同一 平面 中 得 到 配 准 图 像 。 以 图 像 所 在 的 平 
面 为 参考 平面 ,对 于 图 像 j 中 的 任意 一 点 x,y ), 其 在 图 
像 中 的 投 昧 *, y 站 式 (5) 求 得 。 

Gy 了 =H (x,y, D) (5) 
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图 1 提取 sift 特 征 点 "" 
Fig.l Extracting sift features. 


图 2 是 3 张大 小 相同 的 图 像 配 准 结果 。 在 配 准时 ， 
从 右 往 左 先 将 两 张 图 配 准 ,再 将 配 准 结果 与 第 3 张 图 进 
行 配 准 。 从 图 中 可 以 看 出 3 张 图 像 能 正确 配 准 , 只 是 图 
像 出 现 了 不 同 程度 的 形变 ,图 中 虚线 框 表示 3 张 原 图 的 
边界 投影 到 同一 平面 上 的 结果 ,从 图 2B 的 边界 框 可 看 


图 2 未 经 几何 校正 的 图 像 配 准 


Fig.2 Image alignment before distortion correction(from A to OC). 


1.2.2 图 像 配 准 的 预 处 理 Kaynig 等 "针对 TEM 图 像 
畸变 提出 了 用 多 项 式 来 拟 合 非 线 性 校正 场 的 方法 ,并 借 
助 匹配 的 sift 特 征 点 来 估计 多 项 式 参数 。 对 于 任意 点 
(x,y), 其 校正 公式 为 : 
F=G(x, B=, xy ty) 
(0 0， aaa 


其 中 ,d 为 拟 合 校正 场 的 最 高 阶 ,本 文中 4d= 5,/ 是 由 校正 


(6) 


系数 5 构 成 的 向 量 。 
对 于 所 估计 的 校正 场 , 当 所 有 特征 点 经 过 校正 后 ， 


通过 旋转 和 平移 转换 到 同一 平面 坐标 时 ,匹配 的 特征 点 
能 近似 对 齐 。 若 图 像 经 过 仿 射 变换 4 与 图 像 在 同一 
平面 对 齐 , 令 R 和 7 为 4 中 提取 出 的 旋转 矩阵 和 平移 
矩阵 , 则 满足 : 


a 
二 
部 
四 
如 . 
© 


wT Principle Ee 


or ientation 


] 
| 


一 


出 边界 之 间 不 再 完全 平行 或 垂直 ,并 且 从 岁 2C 的 边界 
框 可 看 出 , 随 着 配 准 图 像 数 目的 增加 ,这 种 形变 会 变 
大 。Kaynig 等 “说 明了 形变 的 原因 在 于 聚 光 透 镜 系 统 
对 图 像 产生 的 几何 畸变 。 因 此 ,在 图 像 配 准 之 前 ,应 先 
进行 图 像 畸 变 校正 预 处 理 。 


2 


上 (xyae' | 加 ($xo)pe’ | 


大 
如 3% @ ol 
A | 
(7) 
min by |X®4" XA" | (8) 
| g 


其 中 ,Kk 为 佑 计 校 正 场所 输入 的 图 像 数 ,本 文中 k=9， 
4=0.01。 

图 像 校正 结果 如 图 3 所 示 , 人 A 为 校正 前 的 图 像 ,B 为 
校正 后 的 图 像 ,C 为 本 文 方法 中 提取 的 非 线性 校正 场 
下 。 通 过 图 3C 可 以 观察 到 ,该 校正 方法 仅 使 图 像 的 局 部 

言 息 发 生 微小 的 形变 ,在 校正 后 仍 能 保持 图 像 的 结构 、 
形态 等 细节 信息 。 图 4 为 图 像 经 畸变 校正 后 再 进行 图 
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像 配 准 的 结果 ,将 该 图 中 虚线 所 示 的 边界 线 与 图 2 中 的 
比较 ,可 以 看 出 图 像 经 过 畸变 校正 后 ,能 有 效 减 小 图 像 
形变 。 然 而 ,图 像 进行 校正 及 配 准 后 ,在 图 像 重 伙 区 域 


的 边缘 仍然 存在 明显 的 接 颖 ,如 图 4 中 白色 箭头 所 示 ， 
因而 需 进行 后 处 理 以 消除 接 缝 。 


图 3 图 像 校 正 前 后 的 结果 及 非 线性 校正 场 F 


Fig.3 Results before, after distortion correction and non-linear correction field F. A: Original image; 
B: Image after distortion correction; C: Non-linear correction field F. 


图 4 校正 后 图 像 配 准 的 结果 


Fig.4 Images alignment after distortion correction with white arrows indicating obvious seams. 


1.2.3 图 像 配 准 的 后 处 理 本 文采 用 Poisson 图 像 编 辑 
方法 "将 重 全 区 域 的 图 像 融 合 到 配 准 的 图 像 中 以 去 除 
接 颖 。 如 图 5 所 示 , 令 g 为 重 羡 区 域 图 像 ,v 为 g 的 梯度 
场 ,s 为 两 张 配 准 后 的 图 像 ,0 为 g 融 合 于 s 时 对 应 的 区 
域 ,60 为 Q 的 边界 。f 和 分别 为 定义 在 Q 内 和 Q 外 表示 
灰 度 的 标量 函数 。 

为 确保 消除 接 颖 的 同时 , 重 琶 区 域 的 梯度 在 融合 前 
后 近似 相等 ,构造 以 下 能 量 函 数 : 


mn 外 | ， 款 flo=f |, (9) 
该 方程 的 解 满足 如 下 泊 松 方程 : 
Af=divy, flo=f "|, (10) 


(10) 式 离散 化 的 最 优 解 满足 : 
| 六 有 二 
dsNhnoa geN, {lon qeN, 


其 中 p 为 重建 区 域 f 中 的 任意 点 ,Np 为 p 点 的 四 邻 域 ,q 
为 四邻 域 中 任 一 点 ,6 为 v,,, 在 PG 方向 上 的 投影 
INp| 为 点 p 的 邻 域 点 个 数 。 该 线性 方程 可 以 通过 
Gauss-Seidel 迭 代 方 法 进行 求解 。 

图 6 为 泊 松 图 像 编 辑 方法 消除 图 4 中 接 颖 的 结 
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图 5 泊 松 编辑 方法 示意 图 
Fig.5 Schematic diagram of poisson image editing. f and fr 
are intensity functions inside and outside area respectively. 


从 图 中 白色 第 关 所 指 处 可 以 看 出 该 方法 能 够 很 好 地 去 
除 接 颖 。 
1.3 多 幅 图 像 的 拼接 

若 答 入 n 张 有 序 的 TEM 图像 , 先 提取 相 邻 9 张 图 像 
的 sift 特 征 点 ,进行 特征 点 的 匹配 ,根据 匹配 的 特征 点 对 
图 像 进行 几何 畸变 校正 ,并 用 所 估计 的 校正 场 校正 所 有 
的 图 像 。 然 后 ,对 所 校正 的 图 像 进 行 两 张 图 像 的 配 准 ， 
再 进行 配 准 后 的 无 颖 融合 以 得 到 两 张 图 像 的 拼接 。 最 
后 通过 图 像 两 两 拼接 的 方式 得 到 多 幅 图 像 的 拼接 。 算 
法 流程 如 下 。 
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图 6 泊 松 图 像 编 辑 方法 消除 接 颖 


Fig.6 Removing seams using Poisson image editing with white arrows indicatingoriginal seams 


removed. 
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图 7 n 张 图 像 拼接 流程 图 


Fig.7 Flow chart of multiple image stitching. 


2 结果 

本 文 实验 数据 来 源 于 南方 医科 大 学 电子 显微镜 室 
采集 的 具有 30% 重 辣 区 域 的 肾 小 球 横 切面 图 ,该 组 数据 
为 60 张 5000 倍 放大 倍数 的 TEM 图像。 实验 选择 了 目 
前 比较 稼 用 的 图 像 处 理 软件 中 的 拼接 方法 即 侣 i 训 中 
grid/collection stitching "功能 、photoshop cs4 2 中 
photomersge 功 能 以 及 ICE2 ,与 本 文 方法 进行 比较 。 

图 8 为 4 种 拼接 方法 对 一 个 完整 骨 小 球 截面 图 的 拼 
接 结果 ,A~D 分 别 为 人 i 软件 中 grid/collection stitching 
功能 .photoshop 中 photomerge 功 能 ICE 软件 的 拼接 功 
能 以 及 本 文 方法 的 结果 。 图 像 采集 的 范围 是 样本 中 包 
含 肾 小 球 的 矩形 区 ,由 图 8 可 以 看 出 ,各 种 拼接 方法 的 
结果 均 有 一 定 程度 的 形变 一 锯齿 样 边缘 ,而 比较 四 幅 拼 
接 图 边缘 与 白色 边框 的 距离 可 知 :本文 方法 的 形变 最 
小 ,其 主要 原因 是 图 像 配 准 之 前 进行 过 图 像 校 正 。 

在 图 像 的 配 准 与 颖 融合 的 处 理 中 ,图 像 的 部 分 重 夺 
区 域 往往 需要 兼顾 多 张 图 像 的 配 准 与 融合 ,因此 处 理 好 


图 像 ,B~E 图 为 图 8 白色 边框 区 域 对 应 的 放大 图 ,该 区 
域 的 拼接 需要 兼顾 3 张 图 像 ,F~J] 图 为 A~E 图 白色 边框 
区 域 对 应 的 放大 图 。 从 图 中 可 以 看 出 :B 图 整个 区 域 都 
产生 了 伪 影 且 在 底部 矩形 框 所 示 区 域 有 明显 的 不 对 章 ， 
C 图 在 箭头 所 指 区 域 存在 严重 的 不 对 齐 ,D 图 中 没有 出 
现 图 像 不 对 齐 , 但 会 造成 图 像 细 节 信 息 的 错误 ,如 I 图 中 
椭圆 标示 人 处。 

通过 比较 可 以 得 出 :本 文采 用 的 校正 方法 能 够 减 小 


mage 


Poisson 
image 
editing 


Output 


图 8 四 种 不 同 拼接 结果 
Fig.8 Results of four different stitching methods, while 
white rectangular areas Sstitched juggling 3 adjacent 
imageand enlarged results show in Fig.9 B~E. A: Stitching 
result of fiji; B: Stitching result of photoshop; C: Stitching 
result of ICE; D: Stitching result of proposed method. 


整 幅 图 像 拼接 产后 的 形变 ,而 且 经 校正 后 再 配 准 能 够 得 
到 比 生 iphotoshop 更 准确 的 配 准 结果 ,同时 本 文采 用 
的 图 像 融合 的 方法 能 够 避免 人 ji 方法 中 图 像 融合 后 产生 
的 伪 影 以 及 ICE 方 法 中 图 像 细节 信息 的 错误 。 
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图 9 局 部 拼接 结果 
Fig.9 Stitching results of local area. A:Original image; B: Stitching result of fiji; C:Stitching result of 
photoshop; D: Stitching result of ICE; E:Stitching result of proposed method; F~J: Enlargedregions of 
rectangular areasin A~E. The white oval in image I indicatesartifact. 
3 讨论 


3.1 图 像 配 准 误差 的 评价 

本 实验 随机 选取 6 组 相 邻 的 图 像 ,通过 计算 图 像 配 
准 的 误差 作为 对 本 文 方法 的 评价 。 以 误差 函数 即 公 式 
(12) 作 为 量化 评价 的 指标 ,误差 值 越 小 ,说 明 配 准 方法 
越 有 效 。 
(12) 


图 10 是 由 匹配 的 6 组 图 像 中 特征 点 的 误差 制 成 的 
箱 图 ,每 组 中 箱 的 上 下 两 端 各 表示 特征 点 误差 的 1/4 和 
3/4 的 分 位 数 , 箱 中 的 粗 横 线 表示 每 组 误差 的 平均 值 。 
“+ ”表示 异常 点 。 原 图 中 像素 的 边 长 为 23 nm, 本 文 配 
准 误差 范围 是 62.5~187.5 nm。 
3.2 图 像 的 重 受 面积 和 71 对 平均 配 准 误差 的 影响 

图 像 的 重 全 面积 直接 影响 着 特征 点 的 个 数 以 及 计 
算 的 复杂 程度 。 图 11A 反 映 了 图 像 的 重 又 率 对 配 准 结 
果 及 运行 时 间 的 影响 。 实 验 结果 表明 图 像 面 积 的 重 芭 
率 越 大 ,图 像 平 均 配 准 误差 越 小 但 是 配 准时 间 越 长 ,在 


A 
350 一 一 一 Average registration error 20 
< 300 一 一 一 Time 
Ee 
3 
总 
0 


一 已 IIP 
nn SD uu 
SO OO OO 


un 
SS 


Average registration 
上 
S 
SS 


i . 0 
0 0 0.2 0.3 0.4 
Overlapped area rate 
11 图 像 重 又 面积 对 平均 配 准 误差 及 时 间 的 影响 


Fig.11 Meters on average alignment error. A: Overlapped rate on average alignment error and time; B: Lamda on average 


alignment error. 
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10 6 组 图 像 中 特征 点 的 配 准 误差 (px) 


Fig.10 Alignment errorsbetween matched features in 
six group of images. 


重 益 率 超过 0.2 时 ,图 像 的 平均 配 准 误差 变化 不 大 ,但 是 
配 准 时 间 却 显著 上 升 。 综 合 考虑 配 准 误差 以 及 时 间 , 建 
议 图 像 面积 重 酸 率 取 0.2。 
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公式 (7) 中 的 正则 化 项 是 用 来 惩罚 校正 前 后 偏离 太 

远 的 特征 点 对 ,在 4 取 适 当 值 时 ,图 像 能 在 校正 后 保留 图 
像 细节 信息 ,有 助 于 图 像 正确 配 准 。 图 11B 的 实验 结果 
表明 , 当 1 小 于 0.01 时 ,平均 配 准 误差 在 4px 附近 波动 ， 
当 4 大 于 0.01 时 , 配 准 误差 迅速 上 升 。 本 文 建议 /到 0.01。 

在 多 张 图 像 的 拼接 中 ,通过 图 像 的 两 两 拼接 会 造成 
配 准 记 ee 误差 可 以 先 对 所 有 配 
准 参 数 进行 整体 优化 后 再 进行 图 像 配 准 。 这 将 是 本 文 
今后 的 研究 方向 。 

TEM 图 像 能 够 在 高 分 辩 率 下 真实 准确 地 反应 组 织 
和 器 官 的 内 部 信息 ,在 临床 病理 活检 标本 的 分 析 诊 断 中 
起 着 重要 的 作用 。 然 而 ,TEM 图 像 难 以 兼顾 高 分 状 率 
和 大 视野 的 需求 ,为 此 ,本 文 提出 了 一 种 基于 特征 点 的 
TEM 多 幅 图 像 拼 接 方法 。 该 方法 不 仅 能 够 避免 图 像 产 
nt dd 图 像 
埋 息 。 实 验 结果 表明 ,本 方法 能 够 实现 大 视野 高 分 辩 率 
的 TEM 图 像 ,获得 高 分 辩 率 下 的 完整 肾 小 球 ,为 肾 活 检 
病理 诊断 提供 有 价值 的 信息 ,期 待 将 来 可 以 将 此 技 
术 方 法 推广 至 透射 电镜 其 它 临床 活检 标本 的 超 微 病理 
诊断 中 。 
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